Unusual temperature behavior of the entropy of 
the antiferromagnetic spin state in nuclear matter 
with an effective finite range interaction 
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The unusual temperature behavior of the entropy of the antiferromagnetic (AFM) spin state in 
symmetric nuclear matter with the Gogny D1S interaction, being larger at low temperatures than the 
entropy of nonpolarized matter, is related to the dependence of the entropy on the effective masses 
of nucleons in a spin polarized state. The corresponding conditions for comparing the entropies of 
the AFM and nonpolarized states in terms of the effective masses are formulated, including the low 
and high temperature limits. It is shown that the unexpected temperature behavior of the entropy 
of the AFM spin state at low temperatures is caused by the violation of the corresponding low 
temperature criterion. 
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The issue of spontaneous appearance of spin polarized 
states in nuclear matter is a topic of a great current in- 
terest due to its relevance in astrophysics. In particular, 
the scenarios of supernova explosion and cooling of neu- 
tron stars are essentially different, depending on whether 
nuclear matter is spin polarized or not. On the one hand, 
the models with the effective Skyrme and Gogny nucleon- 
nucleon (NN) interaction predict the occurrence of spin 
instability in nuclear matter at densities in the range from 
go to 6pn for different parametrizations of the NN poten- 
tial 

Q-El (ft) = 0.16 mi" 3 is the nuclear saturation 
density). On the other hand, for the models with the 
realistic NN interaction, the ferromagnetic phase transi- 
tion seems to be suppressed up to densities well above 

go M-M- 

Here we continue the research of spin polarizability of 
nuclear matter with the use of an effective NN interac- 
tion. As was shown in Ref. [l3| . in symmetric nuclear 
matter with the Gogny D1S effective interaction the an- 
tiferromagnetic (AFM) spin ordering with the oppositely 
directed neutron and proton spins sets in, beginning from 
some critical density (at zero temperature, g c w 3.8po)- 
At finite temperature, the entropy of the AFM spin state 
demonstrates the unusual behavior being larger than the 
entropy of the nonpolarized state at low enough tem- 
peratures. The main goal of this work is to clarify this 
unexpected moment, utilizing the approximation of the 
effective mass in the single particle spectrum of nucle- 
ons. As will be shown later, this approximation is quite 
successful in reproducing the entropy of a spin polarized 
state for all relevant temperatures. The use of low and 
high temperature expressions for the entropy will allow 
us to get the corresponding conditions in terms of the 
effective masses for comparing the entropies of spin po- 
larized and nonpolarized states. 

Now we stop on the basic equations of the the- 
ory [ll], [l2| ■ Given the possibility of phase transitions to 
the states with parallel and antiparallcl ordering of neu- 
tron and proton spins, the distribution function / can 
be expanded in the Pauli matrices cr,; and in spin and 



isospm spaces 



/(p) = /oo(p)^oTo + /3o(p)ct 3 t (1) 

+ /03(p)o-0T 3 + /33(P)0-3T 3 . 

Expressions for the distribution functions /oo , /30 > /03 1 /33 
read p], Q3 
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Here n{u>) — {exp(w/T) + 1} 1 and 



WnT = £0 + £00 + £30 + £03 + £33 — M"T i 
w pT = £ o + £00 + £30 — £~03 — £33 — MpT ' 

^nl = £0 + £00 — £30 + £~03 — £33 — MraJ.! 

u pI — £0 + £00 — £30 — £03 + £33 — Mpi I 



(3) 



are the branches of the quasiparticle spectrum corre- 
sponding to neutrons and protons with spin up and spin 
down, and their respective chemical potentials 

(t = n,p; (7 =|, {). Under derivation of Eqs. @, ([3]), it 
is assumed that the populations of neutrons and pro- 
tons with spin up and spin down are held fixed. In 
Eq. ([3]), £o(p) is the free single particle spectrum, and 
£00, £30, £03, £ 33 are the Fermi liquid (FL) corrections to 
the free single particle spectrum, related to the FL am- 
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plitudes C/o(k), 1/3 (k) by formulas 

£oo(p) = ^ ^oW/oofa), k 



£30 (p) 

^03(P) 
£33(P) 



(4) 



±^C/ 2 (k)/ 03 (q), 
q 

^-^(7 3 (k)/ 33 (q). 



The distribution functions /oo j /03 , /30 1 /33 , m turn, 
should satisfy the normalization conditions for the to- 
tal density g n + g p = g, excess of neutrons over protons 
Qn ~ Qp = a Q-, ferromagnetic (FM) g^ — q± = Ag-ff and 
antiferromagnetic (AFM) (gnt + Ppi) - (£ , ?4+£?pt) = ^-QU 
spin order parameters, respectively (a being the isospin 
asymmetry parameter, g] = g n \ + g p ] and g[ = g n [ + g P i 
with Q n f, g n ± and g p ^, g P i being the neutron and proton 
number densities with spin up and spin down, respec- 
tively). To check the thermodynamic stability of differ- 
ent solutions of the self-consistent equations ©-([1]), it 
is necessary to compare the corresponding free energies 
F = E — TS, where E is the energy functional and the 
entropy reads 



S = ~y~] ^2 {n(uj Tcr ) lnn(u> T(T ) 

p r=n,p cr=T, I 

+ n(u) Ta ) lnn(u Ta )}, n(uj) = l — n(u). 



(5) 



The single particle energies ([3]) have the following gen- 
eral structure 
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WrCT (k)= w ° CT (k) + C/ rCT (k), w ° ff (k) 
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Mrcr, (6) 



where mo is the bare nucleon mass, U T<7 is the single par- 
ticle potential. Its momentum dependence can be char- 
acterized by the effective mass m Ta (k), defined as 



m 
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dk 
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If to use the quadratic approximation for the single par- 
ticle potential 



U T(7 (k) « U Ta (0) + 



1 dU Ta {k) 
2k dk 
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where kp Ta is the Fermi momentum of nucleons in the 
state (r, a) , then the single particle energy can be repre- 
sented in the form 



2 1,2 



w r<7 (k) 



h'k- 
2m Ta 



+ U Ta (0)-fi TC 



m ra {k Frcr ). (8) 



Within this approximation, all thermodynamic quanti- 
ties can be easily calculated, analogously to the case of 
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FIG. 1: The entropy per nucleon, measured from its value in 
the nonpolarized state, for the AFM spin state as a function 
of temperature at different polarizations. 



a free Fermi gas. Note that in order to get the effec- 
tive mass m Tcr , it is necessary to find the explicit single 
particle potential as a result of the solution of the self- 
consistent equations ©-([1]). 

Further we will consider symmetric nuclear matter 
with the Gogny D1S interaction as a potential of NN in- 
teraction. It was shown in Ref. [l3[ that in this case the 
AFM spin ordering can be realized only among the states 
with the collincar spin ordering. In the AFM spin state, 
Qn] = QpiiQni = Qp], and, hence, ^i nT = /i p j.,/i„| = fi p ]. 
Besides, we have only two different branches in the quasi- 
particle spectrum, w n f = ui p i, and u) n ± = u) p ^, and, 
as a consequence, only two different effective masses, 
?7!,„j = m P i, and m n j = m p f. The neutron II„ = e " T e e " 1 

and proton n p = gpT - Bpl spin polarization parameters 



•jp 



for the AFM spin state are opposite in sign and equal to 



n„ 



_n p = ^Ii^n. 

Q 



In Fig. [TJ the difference between the entropies per nu- 
cleon of the AFM and nonpolarized states is shown as a 
function of temperature at different fixed polarizations. 
The value II = 1 corresponds to the totally AFM polar- 
ized nuclear matter. One can see that for low temper- 
atures the entropy of the AFM state is larger than the 
entropy of the normal state. It looks like the AFM state 
at low finite temperatures is less ordered than the non- 
polarized state. Under a further increase of temperature 
the difference between the entropies changes the sign and 
becomes negative, that corresponds to the intuitively ex- 
pected behavior. 

Fig. [5] shows the difference between the free energies 
per nucleon of the AFM and nonpolarized states as a 
function of temperature for the same polarizations as 
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FIG. 2: Same as in Fig.[TJ but for the free energy per nucleon, 
measured from its value in the nonpolarized state. 

in Fig. [TJ In contrast to the difference between the en- 
tropies, the difference between the free energies preserves 
its negative sign for all temperatures below the corre- 
sponding critical temperature, and, hence, the AFM spin 
state is thcrmodynamically preferable as compared to the 
nonpolarized state for the whole corresponding tempera- 
ture domain. Therefore, the unusual temperature behav- 
ior of the entropy of the AFM spin state doesn't lead to 
the instability of the polarized state at low temperatures. 
Note that if to assume the quadratic approximation for 
the dependence of the free energy per nucleon of the po- 
larized state on the spin polarization parameter poT ]. 

then the curve in Fig. [21 corresponding to II = 1, in 
fact, shows the dependence of the spin-isospin symmetry 
parameter 7 on temperature. Its negative value proves 
the stability of the AFM spin state at temperatures below 
the critical temperature, as clarified above. 

To understand the unusual behavior of the entropy of 
the AFM spin state at low temperatures, we utilize the 
quadratic approximation for the single particle spectrum 
of nuclcons, Eq. ([8]). Note that, adopting this approxi- 
mation, we self-consistently determine also the chemical 
potentials to guarantee the fulfillment of the normaliza- 
tion conditions for the distribution functions. The results 
of the numerical determination of the entropy density, 
based on the exact and approximated forms of the sin- 
gle particle energies, are compared in Fig. [31 It is seen 
that the approximation turns out to be quite satisfactory, 
especially in the region of low temperatures, and the dis- 
agreement does not exceed 7% even for the less favorable 
case of totally polarized matter in the entire temperature 
domain under consideration. Therefore, we can provide 



FIG. 3: (Color online) The density of entropy for the AFM 
spin state as a function of temperature, calculated with the 
exact (solid line) and approximated (dashed line) single parti- 
cle energies at the spin polarization parameter II = 0.6 (left) 
and II = 1 (right). 

the low temperature expansion for the entropy using the 
approximation of the effective mass in the single particle 
energies, as described in detail, e.g., in Ref. [23]. Then, 
requiring for the difference between the entropies per nu- 
cleon of the AFM and nonpolarized states to be negative, 
one can get the condition 

D 1 = ^(1 + u)i + ^(1 - n)* - 2 < 0. (9) 

m* to* 

Here to* is the effective mass of a nucleon in nonpo- 
larized nuclear matter at the corresponding temperature 
and density. The low temperature condition §§§ is valid 

h 2 k 2 

until T/£p n<r <C 1, £f„„ = 2m f " g being the Fermi en- 
ergy of neutrons with spin up (a =f) and spin down 
(a =[). The calculations show that at the given den- 
sity (f? = 0.64 fm -3 ) and polarizations, the correspond- 
ing temperature interval extends approximately up to 
T = lOMcV. Besides, under derivation of the condi- 
tion © it is assumed that the effective masses are tem- 
perature independent. Fig. [4] shows the dependence of 
the effective masses m n -[, m n i on temperature in the tem- 
perature domain, where the low temperature expansion 
holds true. It is seen that the effective masses for par- 
tially AFM polarized nuclear matter at the given polar- 
izations are practically independent on temperature, and 
for totally polarized matter the change in the effective 
mass is about 5% in this temperature domain. There- 
fore, the use of Ineq. ([9]) is quite justified for comparing 
the entropies of the AFM and nonpolarized states at low 
temperatures. 

Fig.Oshows the dependence of the l.h.s. D x of Ineq. @ 
on temperature in the low temperature domain under 
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FIG. 4: (Color online) The effective masses of spin up and spin 
down neutrons in units of the bare nucleon mass as functions 
of temperature in the low temperature domain at different 
polarizations. 

consideration at different polarizations. It is seen that 
the quantity D\ is positive and increases with the spin 
polarization. This explains the unexpected behavior of 
the entropy of the AFM spin state, being larger than that 
of the nonpolarized state at low temperatures. 

Note that at higher temperatures the entropy of the 
AFM spin state becomes smaller than the entropy of the 
nonpolarized state. To explain this, we can again use 
the approximation of the effective mass in the single par- 
ticle energies for getting the high temperature expres- 
sion for the entropy. If ewA 3 ^ < 1 {\ n<J = J l^K? 
is the thermal wavelength of neutrons with spin up and 
spin down), then the condition for the difference between 
the entropies per nucleon of the AFM and nonpolarized 
states to be negative is 

Fig. [5] shows the dependence of the quantity on tem- 
perature in the high temperature region at the given 
density and polarizations. One can see that the condi- 
tion (fT0|) is fulfilled, and, hence, the entropy of the AFM 
spin state turns out to be smaller than that of the nonpo- 
larized state. Note that from the derivation procedure of 
the high temperature expression for the entropy it follows 
that the effective masses in Ineq. (fTU| can be temperature 
dependent. 

Thus, using the approximation of the effective mass, 
it is possible to explain both the low temperature and 
high temperature peculiarities of the entropy of spin po- 
larized nuclear matter with the D1S Gogny interaction. 
In Ref. fl5| , it was found that the entropy of spin po- 



FIG. 5: The differences Dt (left) and D 2 (right) in Ineqs. © 
and (|10|) as functions of temperature at different polarizations. 



larized neutron matter with the Skyrme effective inter- 
action above some critical density gs is larger than that 
of nonpolarized matter. As a consequence, the critical 
density for the appearance of a spin polarized state de- 
creases with temperature, contrary to the intuitive sug- 
gestion. However, there is an important difference be- 
tween the cases with the finite range Gogny and zero 
range Skyrme forces: While for the Skyrme interaction 
the effective masses in Eq. (|7|) arc momentum and tem- 
perature independent, for the Gogny interaction they do 
depend on momentum and temperature. By this rea- 
son, the difference between the entropies of polarized 
and nonpolarized states is positive for all temperatures 
at densities above gs for the Skyrme interaction. For 
the Gogny interaction, this difference changes the sign 
from the positive one at low temperatures to the negative 
one at high temperatures. The critical density gs, above 
which the difference of the entropies becomes positive at 
low temperatures, depends on polarization, for example, 
g s (U = 0.3) « 0.07 fm~ 3 , g s (U = 0.6) w 0.08 fa" 3 , and 
ps(II = 1) « 0.17 far 3 . For comparison, for totally po- 
larized neutron matter with the SLy4 Skyrme interaction 
g s ~ 0.15 far 3 I53. Since the entropy of the AFM spin 
state is smaller than that of nonpolarized matter at high 
temperatures, the critical density for the appearance of 
the AFM state increases with temperature [14j , in agree- 
ment with the intuitive considerations and contrary to 
the scenario with the Skyrme interaction (l5| . 

In summary, it has been shown that the entropy of 
the AFM spin state in symmetric nuclear matter with 
the Gogny D1S interaction demonstrates the unusual be- 
havior, being larger at low temperatures than the en- 
tropy of nonpolarized matter. By comparing the free 
energies of polarized and nonpolarized states, it has been 
clarified that this unconventional temperature behavior 
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doesn't lead to the instability of the AFM state. This en- 
tropy peculiarity has been related to the dependence of 
the entropy on the effective masses of nuclcons in a spin 
polarized state, which for the finite range D1S Gogny 
interaction do depend on temperature. The correspond- 
ing conditions for comparing the entropies of the AFM 
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